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INTRODUCTION 


The  program  described  in  this  report  is  intended  as  a design  aid 
to  give  the  engineer  a tool  to  establish  the  statit  loads  and  moments 
on  a sabot  segment  while  under  the  in-bore  environment.  The  physical 
characteristics  of  the  sabot  segments  are  also  computed.  Addition- 
ally, from  the  input  data,  a tape  is  prepared  by  the  CDC  6500  com- 
puter for  drawing  sabot  cross  sections  on  the  CALCOMP  570  digital 
plotter,  together  with  a NANCY  digital  printer  plot. 


The  computation  scheme  is  capable  of  handling  multiple  components 
of  differing  densities,  but  the  forces  and  moments  are  treated  as 
though  the  entire  sabot  is  a single  rigid  body.  The  user  must  des- 
cribe the  cross  section  of  the  segment,  and  predetermine  the  manner 
in  which  the  pressure  is  distributed  over  the  surface  of  the  sabot. 
Also,  the  weight  and  first  moments  are  computed  utilizing  signed 
( + 6 -)  densities  to  describe  the  geometry.  The  geometry  must  be  a 
body  of  revolution. 


GENERAL  INSTRUCTIONS 


This  sabot  program,  written  in  Fortran  IV,  consists  of  two  sepa- 
rate sub-programs,  each  having  its  own  input.  The  first  sub-program, 
occurring  in  the  program  listing,  computes  overturning  force  moments 
on  the  sabot  as  a function  of  propellant  pressure  acting  on  external 
sabot  surfaces,  such  as  the  sabot  base  and  loading  of  rear  support 
and  obturation  ring,  etc.  The  remaining  sub-program  computes  over- 
turning moments  as  a sole  function  of  propellant  pressure  acting 
between  internal  sabot  parting  surfaces. 


Instructions  pertaining  to  each  of  these  sub-programs  are  given 
in  succeeding  paragraphs. 


With  regard  to  the  first  sub-program,  a physical  description  is 
accomplished  by  dividing  the  cross  section  of  the  sabot  into  elements, 
each  element  being  originated  at  an  inflection  point  in  the  geometry. 
Each  component  is  broken  down  separately  into  these  elements,  which 
are  entered  as  data  in  x and  y coordinates  in  an  ordered  sequence 
following  the  contour  of  the  components  in  a clockwise  direction.  1 
The  order  is  significant  to  maintain  a system  to  ascertain  the  direc- 
tion of  the  pressure  force  vector.  It  is  also  necessary  to  intro- 
duce signed  densities  to  distinguish  a "forward  surface"  and  a 
"rear  surface".  A zero  density  "end  point"  for  components  is  also 


‘This  portion  of  the  program  was  written  to  accept  a maximum  of 
"5  sets  of  coordinates. 
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used  to  delineate  separate  components.  An  example  of  a two-component 
sabot  divided  into  elements  is  shown  in  Figure  2,  page  22.  The  data 
are  shown  in  tabular  foim  on  page  3.  Note  that  the  concept  of  the 
negative  density  is  associated  with  a "rear  surface,"  and  that  each 
surface  is  described  in  the  element  number  which  originates  that 
surface. 

The  clockwise  flow  of  data  refers  to  the  overall  direction  around 
the  component  and  not  the  direction  of  a point  tracing  an  individual 
line;  points  17  through  28  are  ordered  by  the  overall  direction  around 
the  sabot  body,  though  the  surface  is  curved,  such  that  a point  trac- 
ing the  curve  progresses  in  a clockwise  direction.  The  final  end 
point  is  entered  separately,  and  should  be  the  same  as  the  first  point 
of  the  final  component  in  order  to  establish  a closed  figure.  This 
separate  entry  serves  the  same  purpose  as  the  artificial  "zero  den- 
sity" points  entered  within  the  body  of  the  data. 

A moment  center  is  also  entered  with  x and  y coordinates.  This 
point  is  chosen  as  the  point  of  restraining  force,  usually  located 
on  the  projectile.  This  moment  center  is  the  "hinge  point"  of  a 
segmented  ring- type  sabot.  Since  this  point  will  offset  the  results 
of  the  moment  balance,  it  must  be  chosen  with  care  in  order  to  ac- 
curately represent  the  actual  conditions.  Also,  it  should  be  noted 
that  this  moment  center  is  a single  point  on  a three-dimensional 
figure.  Care  must  be  exercised  to  allow  for  curvature  in  selection 
of  the  moment  center.  In  some  geometrical  construction  the  rotation 
may  occur  about  a shifting  line  of  multiple  points. 

A radius,  with  the  x and  y coordinates  of  the  radius  center,  is 
included  for  segments  which  are  curved  following  the  inflection  point. 
These  values  are  entered  as  zeros  for  straight  line  segments.  The 
next  data  item  is  the  density,  entered  as  positive  for  a forward  sur- 
face and  negative  for  a rear  surface  (a  rear  surface  faces  the  X = 0 
coordinate).  The  next  two  entries  describe  the  pressure  field  to 
which  the  surface  following  the  inflection  point  is  exposed.  The 
first  of  these  (PV)  is  a factor  bv  which  the  chamber  pressure  (en- 
tered later)  will  be  multiplied  to  gi'  the  pressure  to  which  the 
surface  is  exposed  (usually  1.0  or  O.u,.  The  second  code  (LP) 
chooses  two  options  in  depicting  the  pressure  distribution  on  the 
surface.  A code  of  1 denotes  a constant  pressure  equal  to  the  chamber 
pressure  times  the  pressure  factor.  A code  of  2 will  distribute  the 
pressure  in  a linear  manner  between  the  pressure  given  by  this  entry 
and  the  pressure  given  for  the  next  point. 
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There  must  be  one  card  for  each  inflection  point  in  this  portion 
of  the  deck.  The  next  set  of  two  cards  enters  the  moment  centers 
and  the  end  points  on  one  card  and  the  number  of  segments,  the  chamber 
pressure  (psi),  the  spin  (rev/sec),  a code  to  determine  the  output, 
and  up  to  60  characters  of  notes  or  descriptive  information  to  be 
included  in  the  printout.  An  output  code  of  ”0"  will  print  all  data 
points,  "1"  will  print  the  input  and  pressure  distribution,  "2"  will 
print  only  the  forces  and  moments.  These  last  two  cards  must  appear 
as  a set  for  each  run  to  be  conducted. 

The  following  paragraph  pertains  to  input  instructions  for  the 
final  portion  of  this  program.  A physical  description  of  the  in- 
ternal gas-pressurized  area  (defined  by  obturation  zones)  between 
sabot  parting  surfaces  is  accomplished  in  a manner  similar  to  that 
described  aboie.  The  bounded  planar  area  is  broken  down  into  ele- 
ments, each  element  being  originated  at  an  inflection  point  in  the 
geometry.  Each  element  is  entered  as  data  in  S and  T coordinates  in 
an  ordered  sequence  following  the  boundary  of  the  gas-pressurized 
area.2  The  final  end  point  is  the  same  as  the  first  point  entered  to 
complete  the  closed  figure.  Therefore,  no  separate  entry  is  needed 
for  the  end  point.  There  is  no  moment  center  entry,  as  the  one 
entered  in  X and  Y coordinates  in  the  first  portion  of  this  program 
is  used. 

Data  organization  of  batch  processing  a job  requires  a data  card 
deck  in  a format  similar  to  that  shown  on  Figure  3,  page  23.  Multiple 
runs  are  possible  on  the  same  configuration.  For  instance,  several 
sets  of  ballistic  conditions  may  be  run,  or  the  moment  center  may  be 
varied,  or  both  of  these  may  be  combined.  The  initial  card  initiates 
reading  by  giving  the  number  of  runs  and  the  number  of  inflection 
points  in  the  main  data.  This  is  followed  by  the  data  cards  which 
describe  the  cross  sectional  geometry  by  giving  the  X and  Y coordi- 
nates of  each  inflection  point.  This,  in  turn,  is  followed  by  data 
cards  which  describe  the  planar  area  geometry  by  giving  the  S and  T 
coordinates  of  each  inflection  point.  Finally,  a set  of  two  cards 
is  provided  for  each  of  the  respective  program  runs.  These  cards 
give  moment  centers  and  descriptive  text.  An  example  of  data  cards 
giving  four  program  runs  is  shown  in  Sample  C3se  Figure  3 on  page  23. 


2This  portion  of  the  program  was  written  to  accept  a maximum  of 
20  sets  of  coordinates. 
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SCOPE  CONTROL  CARDS  USED  IN  THIS  PROGRAM 


The-  to’ lowing  control  cards  are  used  to  call  NANCY  print  plot 
and  CALC  '’  plot  subroutines  from  permanent  files  stored  in  Picatinny 
A renal’s  CDC  6600  system  (Ref  1,  2 and  3): 

J Ob  , • » • 


Request,  Tape  77,  NT,  S.  Plotape  - your  name 
Attach  ,NAN.  .NANCY,  CY=2,  SD=16,  MR=1,  ID=RANDERS) 

’ oad  (NAN) 

LHC  (Your  program) 

is  a Fortra  : deck  which  writes  a printer  plot  on  Tape  6 and 
a ■ "?MP  plot  on  Tape  77.  The  plot  routines  (Ref  1)  also  require 
that  OUTPUT  be  included  on  the  program  card.  Memory  required  to  load 
is  about  20,000  (octal)  CM  cells. 


INPUT  DATA  CARDS 

1st  Card:  Header  card  (only  fixed  point  entries  on  this  card) 

Column  1,  number  of  computer  runs 

Columns  2 and  3,  number  of  "inflection  points"  about 
cross  section 

Columns  4 and  5,  number  of  "inflection  points" 

about  pressurized  plan;  area  of 
sabot  parting  planes. 

2nd  Group:  Elements  of  sabot  cress  section  (floating  point  except 

for  Column  40  fixed  point  entry) 

Columns  1-6,  X coordinate  of  cross  section 

Columns  7-12,  Y coordinate  of  cross  section 

Columns  13-18,  length  of  radius  of  curvature 


j 


l 


i 
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Columns  19-74,  X coordinate  of  radius  of  curvature 
origin 

Columns  25-30,  Y coordinate  of  radius  of  curvature 
origin 

Columns  31-36,  density 
Columns  37  39,  pressure 
Column  40,  pressure  code 

3rd  Group:  Elements  of  pressurized  planar  surface  on  sabot 

part.ug  planes  (floating  point  entries) 

Columns  1-6,  S coordinate  of  planar  surface  element 

Columns  7-i2,  T coordinate  of  planar  surface  element 

Columns  13-18,  radius  of  curvature  magnitude 

Columns  19-24,  S coordinate  of  radius  of  curvature 
origin 

Columns  25-30,  T coordinate  of  radius  of  curvature 
origin 

4th  Card:  Sabot  segment  moment  center  and  end  points  for  cross 

section  (floating  point  entries) 

Columns  1-6,  X coordinate  of  moment  center 

Columns  7-12,  Y coordinate  of  moment  center 

Columns  13-18,  X coordinate  of  end  point 

Columns  19-24,  Y coordinate  of  end  point 

5th  Card:  Sabot  data  input  and  data  callout  (floating  point 

except  for  Column  21  fixed  point  entry) 

Columns  1-2,  number  of  sabot  segments 

Columns  3-8,  propellant  pressure  (psi) 

Columns  9-12,  diameter  of  -saboted  round  (mm) 


I 


Columns  13-16,  weight  of  subprojectile  (lb) 

Columns  17-20,  spin  of  round  (rev/sec) 

Column  21,  If  "0"  prints  all  data  points 

If  "1"  prints  input  and  pressure  distribution 

If  "2"  prints  only  forces  and  moments 

Columns  22-71,  Data  statement  such  as,  "sabot  test  for 
moment  about  front  edge" 


EQUATIONS 


The  following  equations  are  used  in  this  program: 

The  area  of  the  trapezoidal  element  a,  b,  c,  d (Fig  1)  is: 

AREA  = [2Y(n)  sin  (DBM  2Y(n+1)  sin  (D6)l  x [|Y(n)  - Y(n+1)|Z  + 

|XCn)  - Vi/] 

The  centroid  location  in  the  x-direction  of  the  trapezoidal 
element  (Fig  1)  is: 

XCENT  = [X(n)  - X(n+1)]x  [2Y(n)  sin  (DB)  * 4Y(n+J)  sin  (DB)]/ 

3I2Y(n)  sin  <D8)  * 2Y(n>J)  sin  CDS})  * X(ntl) 

for  case  of  Y,  > Y, 

(n)  (n+l) 

The  weight  f the  prism  element  formed  by  the  "AREA"  multiplied 
by  the  rectangle  element  e,  f,  g,  h (Fig  1)  is 

^ y • y n 

WT  = p x AREA  x cos  j arc  tan  |-&H x (XCENT) 

X(n)  ‘ X(n+1)  J 

p is  density 


The  pressure  forces  acting  on  the  trapezoidal  element  (Fig  1), 
along  X and  Y directions,  are 

F = AREA  x PS  X cos  0 
x 

Fy  = AREA  x PS  x sin  0 
PS  is  pressure  vector. 


IDENTIFICATION  OF  SYMBOLS 
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BETAR 

N 

NJAY 

NN 


NNJAY 


NRUN 

P(K) 

XCENT 

X(K,J) 

S(K,J) 

THETA(K.J) 

T(K,J) 

YCENT 

Y(K,J) 

YVAR 


Angle  in  radians  between  parting  sabot  planes 
Number  of  inflection  points  about  sabot  cross  section 
End  inflection  point  around  sabot  cross  section  to 
close  loop 

Number  of  inflection  points  about  pressurized  zone 
of  sabot  parting  surface 
End  inflection  point  around  pressurized  zone  of 
sabot  planar  parting  surface 
Number  of  computer  runs 
Pressure  factor 
Centroid  in  X-direction 

X coordinates  for  inflection  points  about  sabot 
cross  section 

S coordinate  for  inflection  points  about  pressurized 
zone  on  parting  planes 

Angle  between  the  horizontal  and  any  pressure  vector 
T coordinates  for  inflection  points  about  pressurized 
zone  on  parting  planes 
Centroid  in  Y direction 

Y coordinates  for  inflection  points  about  sabot 
cross  section 

Dummy  integration  variable 
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Fig  1 Geometry  of  Sabot  Cross  Section 
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APPENDIX 


SAMPLE  CASE 


The  following  example  should  give  the  user  a good  idea  of  input 
format  and  the  usual  output  from  this  program. 

This  sample  case  is  illustrated  by  a drawing  of  the  cabot  sec- 
tion (Fig  2)  with  points  of  inflection  shown.  This  is  followed  by 
coding  sheets  (Fig  3)  for  its  input  and  a listing  of  all  its  output 
(Fig  4).  Finally,  the  NANCY  printer  plot  (Fig  5)  and  CALCOMP  digital 
plot  (Fig  6)  are  included  as  they  would  appear  as  part  of  the  com- 
plete data  printout. 


Fig  3 Input  data  cards  format 
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Fig  3 (Continued) 


Fig  3 (Continued) 


Fig  4 (Continued) 
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Fig  4 (Continued) 


Fig  4 (Continued) 
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Fig  4 (Continued) 


Fin  5 NANCY  Drinter  plot  of  sabot  cross  section 


Fig  6 CALCOMP  digital  plot  of  sabot  cross  section 


